Magnetic field dependence of tunnel couplings in carbon nanotube quantum dots 
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By means of sequential and cotunneling spectroscopy, we study the tunnel couplings between 
metallic leads and individual levels in a carbon nanotube quantum dot. The levels are ordered in 
shells consisting of two doublets with strong- and weak tunnel couplings, leading to gate-dependent 
level renormalization. By comparison to a one- and two-shell model this is shown to be a consequence 
of disorder-induced valley mixing in the nanotube. Moreover, a parallel magnetic field is shown to 
reduce this mixing and thus suppress the effects of tunnel-renormalization. 
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Confined states and their tunnel couplings to the leads 
are the basis for a wealth of intriguing phenomena ob- 
served in quantum dot nanostructures. In carbon nan- 
otubes the states appear in a particularly simple arrange- 
ment of near four-fold degenerate shells stemming from 
spin and orbital degrees of freedom [l','^^. The splitting of 
each quartet into two doublets is well understood within 
a single particle model including disorder-induced valley 
mixing and spin-orbit interaction [sl, 0] • Less is known 
about the tunnel couplings, but experiments show that 
they can be different for the two doublets within a shell 
p. In a clean nanotube, such behavior would be incon- 
sistent with time-reversal symmetry. 

In this Letter, we show that the relevant doublets in- 
deed may have different lead couplings due to orbital 
disorder, and that this tunnel coupling asymmetry is 
tunable by a parallel magnetic field B\\. At 5|| = OT, 
strong and weak Kondo effects in the charge stability di- 
agrams and gate-dependent inelastic cotunneling thresh- 
olds indicate different doublet couplings [5] . At finite B\ \ , 
Coulomb peaks exhibit B\\ dependent level broadenings 
indicating field tunable tunnel couplings, which is further 
corroborated by the disappearance of gate-dependent in- 
elastic cotunneling lines for equal couplings. A single 
shell model including disorder, expressed as an intra-shell 
coupling Akk' between the orbital states K and , is 
shown to account well for the observed coupling behav- 
ior at low fields. At larger the coupling strengths of 
the doublets even interchange and show non-monotonic 
behavior, which can be explained by a two-shell model 
including an inter-shell orbital coupling A kk' 12 [6]. By 
applying the functional renormalization group (fRG) to 
the inner four levels of the two-shell model, it is shown 
that the linear conductance versus gate voltage Vg^te is 
in good agreement with the measurement for all B\\ 
including correlation effects such as Kondo effect. 

Figure [TJa) shows a bias spectroscopy plot obtained 



by measuring the differential conductance at tempera- 
ture T = 140 mK of a 400 nm long AuPd contacted single 
wall carbon nanotube device versus source-drain bias and 
backgate voltage (see inset in Fig. [TJa)) [4]. In the V^ate 
range shown, quantum states belonging to two near four- 
fold degenerate shells in the valence band are filled. This 
leads to the characteristic pattern [1] of three smaller 
faintly visible diamonds followed by a big (truncated) di- 
amond as Vgg,te is increased. Inside each diamond the 
nanotube holds an integer number of electrons, indi- 
cated by the additional electron number for filled shells. 
For odd occupancy of the shells, zero-bias conductance 
ridges are observed, which is a manifestation of the well- 
known spin-half Kondo effect [8, 9]. Interestingly, the 
Kondo ridges are broad and narrow for + 1 and N -\-3 
electrons in the shell, respectively. A gate-dependent step 
in conductance at small finite bias (especially pronounced 
for filling N -\- 2 and A/" + 3 in the small diamonds) re- 
veals the splitting of the four states within a shell into 
two doublets via inelastic cotunnehng; i.e., at bias volt- 
ages corresponding to the doublet energy difference, an 
electron may tunnel into the excited doublet provided 
that an electron tunnels out from the ground state dou- 
blet Both features can be understood by assum- 
ing that the two doublets have different couplings to the 
leads: a strongly (weakly) coupled ground (excited) state 
results in a strong (weak) Kondo resonance, while gate- 
dependent level renormalizations of the doublets are due 
to charge fluctuations between the dot and the leads (see 
SM) [5]. Until now it has remained unclear how the two 
doublets could be differently coupled. 

To understand the doublet states in more detail, the 
zero bias conductance of a single shell is measured as a 
function of B\\ and results are shown in Fig. [Ifb) (gate 
range in Fig. [TJa) marked by horizontal dashed line). 
The upwards and downwards energy shifts of pairs of 
Coulomb peaks with increasing field are consistent with 
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FIG. 1. (color online), (a) Charge stability diagram of two nanotube shells showing a pattern of strong/weak Kondo effect for 
odd occupation and gate-dependent inelastic cotunneling thresholds. Inset: backgated two terminal carbon nanotube device, 
(b) Linear conductance versus B\\ for the Vgate region marked by the dashed yellow line in (a), (c) Shell energy diagram showing 
the eigenenergies in the presence of an orbital coupling A^k'- (d,e) Level energies and level-lead couplings of a simple spin-less 
model. Due to Akk' (with orbital phase 00 = tt), the eigenstates are now differently coupled to the leads at 5|| = 0, while 
becoming more equally coupled as B\\ is increased, (f, g) Simple picture showing that the electron probability distributions for 
the two eigenstates are different /similar at zero/large field. 



expectations from a simple nanotube shell model includ- 
ing intra-shell orbital coupling Akk' [UJ with energies as 
shown in Fig.[T]^c) (see SM). The lower (red/orange) 
and upper (blue/cyan) lines are the resulting eigenstates 
corresponding to a finite intra-shell coupling. Moreover, 
the widths of the Coulomb peaks (increasing with tun- 
nel coupling) emerging from the broad Kondo ridge be- 
come narrower as the field is increased (diamonds), while 
a broadening is observed in case of the narrow Kondo 
ridge (squares) [l2[. The tunnel couplings to the leads 
are thus observed to be B\\ dependent. 

The observed strong and weak tunnel coupHng behav- 
ior of the doublets is a consequence of the intra-shell 
orbital coupling as can readily be understood from a 
simple spinless model with level energies shown in Fig. 
WiA) (see SM). Time-reversal symmetry demands that 
I^kI = I^K'I = ^5 but in the presence of a finite intra-shell 
orbital coupling, Akk' = \^KK'\e'^^° ^ so that the new 
eigenstates \b/a) ^ ^e*^o|i^) + at 5|| = are tun- 
nel coupled by t(l ^ 6*"^°), respectively (see SM) [HI. For 
(j)o = TT and 5|| = 0, one of these states completely de- 
couples from the leads, while the original eigenstates are 
restored and couplings become equal for B\\ ^ \^kk'\ 
(see Fig. [U^e)). With this picture in mind, the decrease 
(increase) in coupling versus field of the strongly (weakly) 
coupled doublets is readily understood. The model gives 
a microscopic picture of why the tunnel couplings of the 
two doublets can be different. It does not, however, ad- 



dress the issue of conserved quantum numbers required 
for the SU(4) or orbital Kondo effect SHQ. Figure 
[U^f) shows that, the electron density (and therefore tun- 
neling amplitude) of the two doublets may be different 
depending on the exact tunneling site on the circum- 
ference. As the field is increased, the respective electron 
densities become similar (pattern changes from standing 
waves to plane waves), consistent with equal tunnel cou- 
plings of the doublets (see Fig. [TJg) and SM). 

A striking consequence is the tunability of the level 
renormalization by B\\. Figure [2fa) shows a zoom in 
the stability diagram of the left shell in Fig. [TJa) at 
B\\ = where the onset of inelastic cotunneling processes 
is clearly observed (see the arrow). The gate dependence 
of the onset has been explained by tunneling renormal- 
ization due to a strongly and a weakly coupled doublet 
[H, [m . Knowing that the intra-shell orbital coupling may 
induce differently coupled doublets, we can now further 
test both this simple model and the many-body origin of 
the gate dependence. As B\\ is increased, the couplings 
to the doublets become equal, and as a result the inelas- 
tic cotunneling threshold is expected to gradually become 
gate independent. This is indeed observed in the stability 
diagrams at finite where the gate dependence dimin- 
ishes for B\\=1T (Fig. 2(b)) and is absent for B\\ = 2T 
(Fig. 2(c)). The similar Coulomb peak widths in the 
linear conductance trace at 5|| = 2T indicate that the 
coupHngs within the quartet become nearly identical at 
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FIG. 2. (color online), (a-c) Charge stability diagrams at 
B\\ = 0, 1, 2T showing that the gate dependence (originating 
from differently coupled doublets) of the inelastic cotunneling 
threshold vanishes as the parallel field is increased in accor- 
dance with the field dependence of the doublet coupling asym- 
metry (see (d)). (d) Extracted lead couplings from Coulomb 
peak widths (Fave at B| | = for the lower doublet is a rough 
estimate, since the fitting formula is not valid in the Kondo 
regime - see SM). 



this field. The B\\ dependence of the doublet couplings, 
estimated from the average width of the two Coulomb 
peaks belonging to each doublet, is clearly revealed in 
Fig. M^d) [17] (see SM). 

A more careful analysis of the Coulomb peak widths 
reveals that the single-shell model does not fully account 
for all features. At large fields (> 3T), the coupling of 
the weakly coupled doublet becomes (even) stronger than 
that of the strongly coupled doublet (Fig. 1(b)), while 
the model predicts equal coupling (see Fig. [D^e)). To 
fully understand the experimental behavior of the tunnel 
couplings at large fields, a two-shell model taking into ac- 
count couplings between levels stemming from different 
shells must be invoked [l8|. Figure [H^c) shows the linear 
conductance versus B\\ and Vg^te for the gate range shown 
in Fig. [IJa). An almost identical behavior of Coulomb 
peaks belonging to the two shells is seen, indicating sim- 
ilar tunnel coupling doublet asymmetries [12]. Further- 
more, at larger fields, kinks in the Coulomb peaks and 
slanted conductance ridges in the Coulomb blockade for 
filled shells (e.g. filling 4) indicate level crossing behavior 
giving rise to the Kondo effect. 

For this strongly coupled device, the level structure is 
best probed by inelastic cotunneling spectroscopy. Fig- 
ures [Hb,e,f) show bias spectroscopy plots versus B\\ for 
different fillings of the shells along the white dashed lines 
in Fig. life), i.e., one hole (3 electrons), a filled shell and 
one electron in the next shell. The observed cotunnel- 
ing features are understood by considering two nanotube 



shells separated by a level spacing AE = E2— Ei^ where 
in addition both shells are split into two doublets due to 
(primarily) intra-shell orbital interaction Akk'i^ Akk'2 
(Fig. ma), see also SM). As B\\ is applied, the doublets 
within each shell split up, resulting in level crossings be- 
tween different shells at elevated fields (see Fig. [3fa)). 
The thin black lines in Figs. [Sfb), O^c), and [H^f) repre- 
senting cotunneling excitation, which were obtained from 
the level spectrum of Fig. [Ifa), well reproduce the ob- 
served features [19]. The cotunneling data also reveal 
the presence of an inter-shell orbital interaction Akk' 12^ 
which couples orbital states with same spin from differ- 
ent shells, and appears as anticrossings at finite field (e.g. 
arrows in Figs. [3fb) and[3l^f)). On the other hand, cross- 
ings for states of different spins (e.g. arrow in [Iff)) indi- 
cate that spin- flip scattering is suppressed [20j. 

The orbital couplings not only modify the internal level 
structure, but also induce B\\ dependent tunnel couplings 
as shown in Fig. [H^d). These effective F's are calculated 
from the eigenstates of the two-shell model using the bare 
tunnel amplitudes. At zero field the doublets arrange in 
strongly and weakly coupled doublets, but as the field 
is increased, the couplings non-monotonically approach 
the original couplings of the shells (see horizontal arrows 
- the valence band shell at larger negative V^ate is more 
strongly coupled as often seen in nanotubes). The non- 
monotonic coupling behavior is also observed in the data 
of Fig. [3l^c). For instance the left most Coulomb peak in 
the red square initially broadens and then becomes very 
narrow at high fields in agreement with Fig. [Sfd) (solid 
blue line). It is also seen that the tunnel couplings of the 
inner two doublets in Fig. [3l^a) interchange at low fields 
(see dotted vertical arrows in Fig.[3l^d)) qualitatively con- 
sistent with Fig. [3l^c), where the widths of the narrowest 
Coulomb peaks (two leftmost in the red square) become 
the broadest as the field is increased. Whereas the quan- 
titative details of Fig. [3fd) depend on the phases of the 
orbital couplings (see SM), the agreement with the exper- 
iment is better assessed by comparing the overall features 
to an fRG calculation based on the level structure and 
tunnel couplings. 

The fRG calculated linear conductance (see SM) [7] in- 
cluding inter /intra- shell orbital couplings, spin-orbit cou- 
plings, and the charging energy Uc is shown in Fig. Wig) 
(compare to square in Fig. [3fc)). Only the relevant and 
correctly modeled inner four levels are kept. After choos- 
ing appropriate intra-shell coupling phases the calcula- 
tion reproduces the weak and strong Kondo resonances, 
the non-monotonic magnetic-field dependent Coulomb 
peak widths, as well as Kondo ridges in Coulomb block- 
ade at finite fields. Interestingly, the finite field Kondo 
ridges are seen to be gate dependent, thus giving rise 



to a Igate controlled spin ground state transition [2l|, |22| 



consistent with different tunnel couplings of the states in- 
volved (see SM) [23]. Note that in the T = calculation 
left-right asymmetric level-lead couplings (corresponding 
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FIG. 3. (color online), (a) Energy levels of a two-shell model with level crossings (see vertical lines) resulting in finite B\\ Kondo 
ridges in (c) for levels having different spins. In contrast, avoided crossings appear at zero/finite B\\ due to intra/inter-shell 
orbital couplings for levels having same spin, (b) Inelastic cotunneling spectroscopy probing level differences for 1 hole. Lines 
are fits (to peaks/dips in d^I/dV^d) based on the two-shell model shown in (a), (c) Linear conductance versus B\\ and Vgate 
allowing the study of level crossings of states belonging to different shells, (d) Orbital couplings give rise to effective tunnel 
couplings of the eigenstates in (a) (horizontal arrows indicate bare shell coupling strengths). Note that we show the sum 
of the couplings to the left and right leads. The small shift of the couplings within a doublet along B\\ is due to spin-orbit 
interaction, (e-f) As in (b), but for a filled shell and 1 electron, (g) fRG calculation based on the inner four levels of (a) using 
the corresponding tunnel couplings of (d) 15]. 



to a higher linear conductance than in the experiment) 
were chosen to anticipate a suppression at finite T [24] . 

In conclusion, we have shown that the doublets formed 
in a nanotube shell in presence of disorder-induced valley 
mixing may have different tunnel couplings to the leads. 
Furthermore, this difference is modified by applying a 
parallel magnetic field. The linear conductance £RG re- 
sults for a two-shell model show good agreement with 
experiments. 
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